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We investigate deformations of Lifshitz holography in (n + 1) dimensional spacetime. After discussing 
I the situation for general Lifshitz scaling symmetry parameter z, we consider z = n — \ and the associated 
^ marginally relevant operators. These operators are dynamically generated by a momentum scale A ~ and 
correspond to slightly deformed Lifshitz spacetimes via a holographic picture. We obtain renormalization 
group flow at finite temperature from UV Lifshitz to IR AdS, and evaluate how physical quantities such 
^ as the free energy density and the energy density depend on log(A^/T) in the quantum critical regime as 
A7T ^ 0. 
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1. Introduction 



One of the most innovative ideas in theoretical physics in recent years is the anti-de Sitter spacee- 
time/Conformal Field Theory (AdS/CFT) correspondence [1], which posits an isomorphism between sym- 
metries of a spacetime and those of matter fields - for instance, the conformal group SO{n, 2) of a 
n-dimensional CFT aries as the group of isometrics of AdSn+i- This idea has been extended to a duality 
between gauge theories and gravitational theories in one larger dimension, yielding a new way to under- 
stand physics. The most exciting aspect of this picture is to provide a technique for obtaining a weakly 
coupled and calculable dual description of strongly coupled matter fields in terms of gravity. 

Several years prior to the conception of the AdS/CFT correspondence in high energy physics, inves- 
tigations on phase transitions of modern materials indicated that a precarious point exists between two 
stable phases of matter such as superconductors and ferroelectrics or ferromagnets in which the temper- 
ature of a system has been driven to absolute zero by the application of some external parameter such 
as pressure or an applied magnetic field. Unlike classical critical points, where critical fluctuations are 
limited to a narrow region around the phase transition, at such 'quantum critical points', the critical 
fluctuations are quantum mechanical in nature and exhibit a generalized scale invariance in both time 
and space. Understanding this puzzling behavior has become a major research effort in condensed matter 
physics. Such systems exhibit universally distinct characteristics upon fanning out to finite temperatures, 
and so the effect of quantum criticality is felt without ever reaching absolute zero. A recent attempt to 
understand quantum critical theory involves extending the AdS/CFT correspondence [2], [3] to condensed 
matter systems, yielding a considerably broader range of scope for gauge/gravity duality. 

The signature scaling property underlying quantum critical theory in 2 + 1 dimensions is 

t^XH, x^Xx (1.1) 

where z is the dynamical critical exponent; z = 1 corresponds to conformal invariance, whereas z ^ 1 
implies an anisotropic scaling invariance. Recently a form of gauge-gravity duality was proposed for z ^ 1 
[4], in which the geometrical dual is obtained from the (asymptotic) metric 

=^ [-^ + ^ + -2 J (1-2) 

which is called Lifshitz spacetime and obviously satisfies 

t XH, r ^ Xr, f Xx. (1.3) 

When z = 1 the metric (1.2) is that of (asymptotic) AdS spacetime and (1.1) recovers the conformal 
symmetry of the CFT. When z = 2, (1.1) restores the scaling symmetry of quantum critical theories. The 
more general anisotropic scaling symmetry, z ^ 1, submerged in the gravity theory and the field theory, 
is the foundation for a Lifshitz spacetime/Quantum Critical Theory (Lifshitz/QCT) correspondence. 

One issue associated with this approach is how to obtain non-trivial spacetimes that asymptote to the 
anisotropic metric (1.2). Two approaches have been considered to this end. First, it is obvious from the 
Einstein equations that an anisotropic energy-momentum tensor could support an anisotropic geometry; 
for example a massive vector field with the appropriate asymptotic behaviour can suffice. An alternate 
approach involves adding higher curvature terms into the Einstein action [5] ; by appropriately tuning the 
different gravitational constants, metrics asymptotic to (1.2) can be obtained. In this paper, we follow the 
first approach, investigating the Einstein action coupled to a massive vector field in (n + 1) dimensions. 

Related to Lifshitz field theory, an interesting feature attracting much recent attention is associated 
with renormalization group flow. The action for the Lifshitz field is 

SLif = \j dTd^x{{dr<Pf - k(V20)2^ , (1.4) 
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where it is clear that the anisotropic scahng invariance has z = 2. By perturbing this action with a 
term — (V(/))^, Lorentz invariance is recovered in a nontrivial way via renormahzation group flow, with 
the perturbation playing the role of a relevant operator. 

In this paper, we consider (n+l)-dimensional Lifshitz spacetime and ((n— l)+l)-dimensional Quantum 
Critical Theory(QCT), and study their holographic duality. While QCT is well described in a 2 + 
1 dimensional context, more general theories of physics including the standard model and gravity are 
implemented in a higher-dimensional context. The success of the AdS/CFT correspondence therefore 
provides motivation to understand the extent to which the broader notions of Lifshitz/QCT duality are 
applicable in higher dimensions, and what different behaviour emerges. Motivated by these interests, we 
especially focus on the marginally relevant operators in the QCT extended to higher dimensions, with the 
goal of understanding their behaviour from the perspective of holographic duality, where these operators 
correspond to the deformed Lifshitz spacetime solutions. 

Previous work in this subject has concentrated on the (2+l)-dimensional case [6] . Here we demonstrate 
that renormahzation group flow from Lifshitz spacetime in the UV to AdS spacetime in the IR generalizes 
to any dimensionality in the marginally relevant case yielding deformations of the pure Lifshitz spacetime. 
From a thermodynamic perspective, we find that physical quantities such as the ratios s/T (entropy 
density to temperature), J- /Ts (free energy density to Ts), and 8/Ts (energy density over Ts) exhibit 
progressively weaker dependence on temperature at sub-leading order in log(A^/T) as dimensionality 
increases. We also find that the maximal flux of the vector field near the horizon grows linearly with 
increasing dimension. 

In section 2, the action, equations of motion, and basic setup are introduced, along with an ansatz 
for which all constants are flne-tuned and normalized for both Lifshitz and AdS spacetime. In section 
3, asymptotic solutions consistent with a marginally relevant operator are derived by bringing in a dy- 
namically generated momentum scale A (assumed very small), which deforms Lifshitz spacetime in the 
high energy regime. In section 4, we carry out holographic renormahzation, rendering the action finite 
by constructing proper counterterms. In section 5, we numerically match our asymptotic near-Lifshitz 
solutions with black hole solutions near the horizon. We then describe the renormahzation group flow, 
and compute physical quantities such as the entropy density s, the free energy density and the energy 
density f for n = 3, 4, 5, 6, 7, and 8. 



2. Einstein Gravity with a Massive Vector fields in (n + 1) dimensional spacetime 

The action for gravity in (n + l)-dimensional spacetime coupled to a massive vector field is described by 



S 



-a 



(2.1) 



-[R + 2A]-- 

where Kn+i = \/87rG„+i in which Gn+i is the (n + 1) dimensional gravitational constant, and H = dB 
and is the (n+1) dimensional coupling constant of the vector field. The equations of the motion are 



1 



R 



1 



g^uR - Ac/, 



1 

gv" 



1 



H^pHf - -.gij.uH + 



1_ 

gv"^ 



Bp,By - -g^ivB' 



(2.2) 



and 



V^,H>"' -jB" = (2.3) 

where 7 is the squared mass of the vector field. For the action to yield solutions asymptotic to those 
having the scaling symmetry (1.3), we require the spacetime metric 



ds^ 



dt" 



dr'' 



dx^ + dy^ + 



(2.4) 
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to be a solution to the field equations, where z is arbitrary. Note that in these coordinates r — )• 
corresponds to the boundary of the spacetime. 

The vector potential yielding a stress-energy supporting this metric is given by 

B = ^^dt. (2.5) 

These ansatz and boundary conditions fine-tune the cosmological constant to be 

(z-l)^ + n(z-2) + n^ 
A = , (2.6) 

and the squared mass and the charge of the vector field to be 

(n — l)z n z — 1 , 

7 = V^. e = — (2.7) 

Regardless of the dimensionality of the spacetime, setting z = 1 in (2.4) yields AdSn+i solution 

ds\,s = '^l\-^ + ^ + ;i j (2-8) 

where the vector potential vanishes. As the cosmological constant has been already fixed due to the 
Lifshitz boundary condition we introduce a scaling constant, a, into the AdS metric and adjust its value 
to be 

a = n{n-l) 

(2-l)2+n(z-2) + n2' ^ ' 

Once we fix the cosmological constant (2.6) with space dimension n and dynamical critical exponent z, 
then those values determine the scaling constant for the AdS spacetime metric. 

In order to describe the renormalization group flow which involves breaking the anisotropy of the 
spacetime by running from the UV Lifshitz to the IR AdS, we employ the ansatz 

ds^ = l^(^-f{r)dt^ + ^+p{r){dx^ + dy^ + ---)y (2.10) 

B = -i^h{r)dt (2.11) 
SO for the Lifshitz spacetime 

Lifshitz :/ = 4t P=^: h= (2.12) 
whereas for the AdSn+i spacetime 



AdS:/=p= (3^"_^^ r V'^^ /, = q. (2.13) 

With (2.10) and (2.11), the equations of motion yield three independent non-linear ODEs for {f{r),p{r)^ ^(^)} 

z(4n - 6)/i(r)2 r/'(r) r'^f'{r)'^ {3n - 5^ f {r)p' (r) {n-2fr^p'{rf r^f"{r)_ 

m JW^^fW 2/(r)p(r) 2^0^ JW ~ 
2zh{r)'^ rp'{r) r'^ f {r)p' (r) r'^p'{r)'^ r'^p"{r) ^ 
f(r) p{r) 2f(r)p(r) 2p{r)'^ p{r) ' 
{n—l)zh{r)'^ r'^h'{rY (n — l)r'^ f {r)p' (r) {n — 2){n — l)r'^p' (r)'^ 

W) W) 2/(r)p(r) 4K0^ = °' ^'^'^ ' 
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where X = ('^ ~ 1)^ + (n — 2)z + z^. We shall rewrite the equations of the motion with the new variables 

p(r) = e/" /(r) = e/" h{r) = k{r)/fi^. (2.15) 

These variables have the added benefit of turning the second order differential equations into first order 
and postponing the determination of rescaling ambiguities on f, p, and h. 
For further simplification, we introduce a new variable 

x{r) = ^4x + 4(n - l)zk{rf - 2(n - l)m{r)q{r) - (n - 2)(n - l)q{rf^ ' . (2.16) 

Putting (2.15) and (2.16) into (2.14) gives 

^ (fi — I) 
rx{r) = —2{n — l)zk{r) ^ q(r)x{r), 

rq'{r) = - zfc(r)2 - ^g(r)2 - ^^-L_x(r)2, 

_ X k{r) in-J) x{r) 1 k{r)x{rY 

In terms of the new variables k{r), q{r), and x{r), Lifshitz spacetime is described by 

Lifshitz : q = -2 x = 2Vz - \yfz, k = ~ \ (2.18) 

whereas for AdS spacetime 



/3n — 4 

AdS : q = -2\ , X = 0, = 0. (2.19) 

V n 

As we are interested in the special case for which the energy-momentum tensor of the gravitational field 
and the operator of the vector field become marginally relevant, the above general description is narrowed 
down to the case satisfying z = n — 1. In this case (2.6) becomes 

~ (3n-4)(n-l) 

A = ^ ^. (2.20) 

Henceforth we deal with only this z = n — 1 case. 

3. Asymptotic Behaviour 

We consider the spacetime slightly thermally heated and so slightly deformed from the pure Lifshitz case, 
restricting our considerations to z = n — 1 for which the massive vector field becomes marginal. Under 
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these assumptions, the general form of the solutions near the boundary r — )• is 



1 + 



{z - 1)2 log(rA) 



+ 



{z - l)i-3z + 2{z - ifX) + 2(1 - 3z)log{- log(rA)) 



2z(z- l)4log2(rA) 



+ 



+ (rA)2Mog"(rA)( /3( 1 + 



2{3z - l)log(-log(rA)) 
z{z — 1)2 log(rA) 



+ 



log(rA) 



+ 



(2^2 -4z + l)- 2{z - lY{2z - 1)A + 2(6^2 - 5z + 1) log(- log(rA)) 



+ 



q{r) 



2^^ 



2z(z-l)2(2z-l)log( rA) 
1 z + 2(z - 1)'^A - 2(3z - 1) log(- log(rA)) 



(3.1) 



{z - 1) log(rA) 



rtvf''\o^\rk)[ /?( 1 + 



+ 



+ a 



2z-l 
1 



2z{z - l)3log2(rA) 
-z(4z2 -7z + 2) + 2(2z - l)(3z - 1) log(- log(rA)) 
z{z - l)2(2z - l)log(rA) 



+ 



(2^2 _ 4^ + 1) + 2(z - \f\ - 2(3z - 1) log(- log(rA)) 



log(rA) 
x{r) = 2y/z - lVi( 1 + 
2^2 



2z{z - l)2log2(rA) 



+ 



(3.2) 



{z - 1)2 log(rA) 



+ 



(z - 1)^A + (1 - 3z) log(- log(rA)) 
(z- l)4log2(rA) 



+ 



2z- 1 



(rA)2Mog^(rA) /3 1 + 



-z(4z2 _ 5^ + 1) + 2(6z2 - 5z + 1) log(- log(rA)) 
z{z - l)2(2z - l)log(rA) 



+ 



+ a 



1 _{2z- if + 2(z - ifX - 2{3z - 1) log(- log(rA)) 
log(rA) 2z{z - 1)2 log2(rA) 



+ 



(3.3) 



where A is a momentum scale, generating a marginally relevant mode, whereas A^ is an energy scale 
with n = z + 1 the spatial dimension. As A — t- the solution recovers the pure Lifshitz spacetime. The 
parameters a and /3 describe other modes of the solution, and A is nothing but a gauge choice [6]. We set 
A = for convenience. 

As we are interested in the high energy regime, we expand by introducing an arbitrary scale ^ and 
write 



log(rA) = log(r^) — lo| 



A' 



In the high energy regime where /z ^ A we have 



1 



log(r^) 



^ 1. 



(3.4) 



(3.5) 



Upon expansion, equations (3.1) ~ (3.3) become 



k{r) 



q{r) 



Vz^ 
-211- 



1 + 



(z-l)2log(^) 



3z{z - 1) + 2z{z - 1)2 log(r/i) + 2(3z - 1) log(- log(f )) 

2Z(Z-I)4l0g2(^) 



+ 



(^-l)log(f 



-z + 2z(z-l)2log(r^) + 2(3z-l)log(-log(f)) 



x{r) = 2Vz - l^/z(l + 



(z-l)2log(^) 



2z(z-l)3log2(^) 

z{z - 1)2 log(r^) + {3z - l)log{- log(^)) 
(z-l)Mog2(^) 



+ 



(3.6) 
(3.7) 

(3.8) 
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Using these solutions for k{r), q{r) and x(r), we employ the change of variables (2.15) and (2.16) in reverse 
to obtain the original form of the solutions 

/ (7z - 4) + 2(3z - 1) log(- log(rA)) (232^ - 142^3 + 152^2 _ 572; + 6) 
(rA)2-(-log(rA))^^ V (z - 1)3 log(rA) 4z(z - 1)6 log2(rA) 

(3z-l)2(5z-2)log(-log(rA)) + (3z-l)3log2(-log(rA)) \ 

z(z-l)6log2(rA) ^"') ^ ' 

P^{- log(rA))^ / {5z - 2) + 2{3z - 1) log(- log(rA)) (312^ - 64^3 + 106^^ - 69z + 14) 



^^^^ (rA)2 z(z-l)3log(rA) + 4z2(z _ 1)6 iog2(,.A) 

(3z3 + 26z2-21z + 4)log(-log(rA)) + (3z-l)2(z-3)log2(-log(rA)) \ 
^ z2(z-l)6log2(rA) ^"V ^ ^ 

where and Pq are constants. Furthermore, in the high energy regime the same expansion for eq. (3.9) 
-(3.10) yields 

, 1 , 7z-4 + 2z(z-l)2log(r//) + 2(3z-l)log(log(^)) , \ 



I 5z - 2 + 2z(z - 1)2 log(r/i) + 2(3z - 1) log(log(^)) 



p{p) = ^ 1 ^ \ \,' ' " + • • • (3.12) 

2;(2: - l)3log(^) ' 



upon rescaling the t and x coordinates to 



^^ 7-1317777 TT^- (3-13) 



4. Holographic Renormalization 

In this section we investigate thermodynamic quantities such as free energy density or energy density at 
an asymptotic boundary of the deformed Lifshitz spacetime. We begin with the definition of the free 
energy density 

F = -T\^gZ = TS,{g) (4.1) 

where 5^ and 5* are respectively the Euclidean action and the metric, and Z is the partition function. 
Upon carrying out a variation of the on-shell action, boundary terms arise, and to cancel these out a 
Gibbons-Hawking boundary term is added into the action. After Euclideanization, the action and the 
metric can be explicitly written as 

5,= / d^^^x^g[-^^{R^2K\-\^-U^^^-B^\\^—^ [ d^x^K, (4.2) 

J V^/^n+r 5^ 4 2 y J 

ds,^ = f (^fir)dr^ + 7^ + Pir){dx^ + V + • • • )) , (4-3) 
where e indicates the Euclidean version of the quantities. 
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Calculating the free energy density (the free energy per unit (n — l)-dimensional spatial volume), the 
Einstein-Hilbert action and Gibbon-Hawking term yield 

J^EH = n lim ri/ f{r)p {r)p{r)^ , (4.4) 

J'gh = 2 1^"^ = 2 ^ ( VW)p{r) ^ (4.5) 

where 'jab is the induced metric on the boundary and K^i, is the extrinsic curvature defined as K^^i, = V^n,y 
in which rii, is the normal vector on the boundary surface. The free energy is F = J (P~^xJ-. However for 
the marginally relevant modes both (4.4) and (4.5) are divergent as the boundary (r — )■ 0) is approached. 
We incorporate boundary counterterms [7, 8, 9, 10] into the action as a remedy for this problem. We 
construct these counterterms as a power series in = B'^B^ [6], so as to satisfy covariance at the 
boundary, obtaining 

•^c.T. = ^r^limV^E^.f-^^'-— (4-6) 
= :f^,l^^yW)piry-^j2Cj(Hrf-^-^y (4.7) 

where we have used {B'^ — {z — 1)/^;) instead of B^, since these must vanish for the pure Lifshitz case. 
The coefficients Cj are not constants but rather a series of the logarithmic functions, with at least three 
needed to eliminate divergences. 

The final expression for the free energy density is 

- lim ^Mr)'^( - + ^ + f C, ik(rf - ^i^V V (48) 



2k„+i2 ^ V Pi'r) f{r) ^ ^ 

To obtain the energy density, we use the definition of the boundary stress tensor to the case in which 
additional non- vanishing boundary fields are present [11]. From the boundary stress tensor, we obtain 
the charge via variation of the on-shell action with respect to the boundary fields; this process in our case 
produces 



SS = Y^'^'^^ab + J^SBa. (4.9) 

Here, however we are dealing not with scalar matter fields but with massive vector fields, and so the usual 
charge defined by 

Q = - [ r-^x^Cahr"'' (4.10) 



where ^/a = ^^xx • • • Izz is the spatial volume element, is a boundary Killing fields, and k^, is the unit 
normal vector to the boundary Cauchy surface, is not conserved. The boundary stress tensor r"* must 
therefore be redefined so as to fix the matter fields in the boundary. Employing the vielbein frame defined 
by 

lab = %ieyl ?7 = diag(±l,l,l,---) (4.11) 

where 

e* = eidx" = y/fdr, e^' = ^ dxi. (4.12) 
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We find that the variation of the free energy density retains its original form, but that replaced 
with T"^, where 

6S = VlTlSet + jHBa, (4.13) 

with 

^ab ^ q-a ba^ ^ab ^ ^ab ^ J_j{a-gb)_ 

The energy density is then given by 

£ = + J'Bu (4.15) 

and the pressure is 

V = -^t\. (4.16) 
Computing the distinct components of iS, we find 

nb 2 6S 1 (T^nb T^ab\ 

2lKn+i^j^^ 'y \ gl z ) gl \ gl z 



(4.17) 



and 



= Jl^lim y7Wp(0'^f-:^2;(r) + %)^jc/M0'-^^^y \ (4.18) 
and other component of ^7" become zero. Putting these together into (4.15) yields 

and 

V = -T. (4.20) 
Imposing finiteness of physical quantities of (4.8), (4.18), and (4.19), the coefficients of the counter terms 
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are found to be 



Co =2{2z - 1) - — + 



(4^4 + 2z^ - - 2z + 1) 



Ci =z + 



{2z -l)iz- 1)3 log2(rA) {z - l)5(2z - 1)2 log3(rA) 
;(-log(rA)) 
-l)log3(rA) '"' 
2^3 z(l4z3 _ ^2 _ + 3) + 4z2(2z - l)(3z - 1) log(- log(rA)) 



4z(3z-l)log(-log(rA)) 

(z- l)5(2z- l)log3(rA) ' ^ ' ' 



{2z -l){z- 1)2 log(rA) 2(2z - l)2(z - 1)4 log2(rA) 

1 /z2(34^4 _ 54^3 _^ 72^2 _ 46z + 9) + 8a(2z - l)2(z - 1)^ 



log3(rA)V 2z(2z- l)2(z- 1)6 

(6z^ + 25^3 - 45^2 + 2lz - 3) log(- log(rA)) 2z(3z - 1)^ log2(- log(rA)) 
(2z-l)2(z-l)6 {z-lf{2z-l) 



+ •••, (4.22) 



_ z'^{l-Zz) z2(15z2 - 14z + 3) a z(3z - l)2(5z - 3) log(- log(rA)) 

^ ~4(2z-l)(z-l) ^ 4(2z- l)2(z- l)3log(rA) ^ log2(rA) 4(2z - l)2(z - 1)5 log2(rA) 

(4.23) 

where the first two are infinite series in l/log(rA) tliat include powers of log(— log(rA)) such that the 
order of the log(— log(rA)) terms do not exceed the order of the 1/ log(rA) terms. It is sufficient for C2 to 
retain terms up to second order in l/log(rA). Note that there exists an ambiguity a in these expressions. 
This ambiguity does not affect numerical evaluation of the free energy density and the energy density 
that we shall later compute, though it does affect reflecting the reaction of the system to changes 
in the boundary Proca field. Our counter term construction (4.21) - (4.23) is minimal; additional terms 
such as C^{B'^ ~ ~ l)/-^)^ or Ca{B'^ ~ ~ l)/-2^)^ would also yield solutions. 

Applying (4.21) - (4.23) into (4.8), (4.18) and (4.19), the physical quantities become 

- K„+i2 ^/^(2z - 1) V {2z-l){z-lY ^J' 

/"-I ( (2z3 _ 42;2 + 4^ _ 1) \ 

K„+i2 V^^(2z - 1) V (2z-l)(z-l)4 ^ ^ 

£_ 1 r-2 /2(20z5 + 18z^-22z3-23z2 + 24z-5) 4(z - l)a\ 
^ ~ ~~ ^71 T^IT^n ^ Z \P- (4.26) 



9v AJn+l V 2(z - l)4(2z - 1)3 

Since the pure Lifshitz solution does not depend on a and /3, we expect in this case that 

7" = f = J* = (4.27) 

regardless of the dimension of spacetime. 



5. Finite Temperature 

In this section, we consider the finite temperature theory by expanding the black hole solution near horizon. 
Our goal is to describe the renormalization group (RG) flow under the marginally relevant modes, and 
to predict behaviour of physical quantities such as the free energy density T and the energy density E 
near criticality in order to provide a way of understanding the quantum phase transition from one phase 
to the critical point via the holographic dictionary. For these purposes, the RG flow corresponding to 
the /iQ-dependent horizon flux is explained by ensuring the RG flow in the zero temperature limit [4] 
A^/T —7- 00 with A ~ fixed. Furthermore J- jTs and E jTs are found as functions of log(A^/r) in the 
near-Lifshitz spacetime containing the large flux. 
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5.1 Expansion and Physical quantities near horizon 

We assume a black hole solution having the form of (2.10) and defined by = 0, and expand the 

solution near horizon r = r+. In this expansion the unknown coefficients arising in every order in the 
f{r),g{r) and h{r) are calculated by the equations of the motion (2.15). Regularity requires gu to have 
a double zero at the horizon, with Qxx remaining nonzero. Under these considerations the functions are 
expressed by 

/(,) . a((, _ (: _ _^y^ M.^ + 14. + 7)..B(3.-2).g ^^ (^^') 

/ (3z - 1)^ - 4/i2 / rV (3z - l)z - 4/i2 / rV \ 

, . ^ rrfh f^ '^\\h(^ , ^ ^ z(-9z^ + 10z + 20) + 8/ig(3z-l) V^ r V 

(5.3) 

where the constants /o and po are associated with scaling ambiguities of the coordinates {t, x, y, • • • }, 
which function as the clock and rulers of the system. Upon fixing these the only variable left in the metric 
is /iQ. Different values of /iq correspond to different black holes and so we have a 1-parameter family of 
black hole solutions. 

Some physical information, for example thermodynamic quantities, can be obtained near the horizon 

of the black hole. At r = r+, the temperature T, can be computed by identifying the imaginary time 

coordinate r with period /3 so as to ensure regularity of the metric, and the entropy density, s, obtained 

A 



from the definition of the entropy, S = -in^ — • We find 



^_r+Jld^f{r) 



(5.4) 



27r V 2 (ir2 

^n— 1 

2vr ^p{r+)^, (5.5) 



r=r+ 

n-l 



where S = J s ^x. 

The horizon fiux, of the massive vector field is 



$ = y VhE ■dA = jcj) cT-^x (5.6) 



where ^ 

^ _ l'^-'^g^r+ / p(r) V dh{r) 

i^n+i V V/(r) dr 
is the horizon fiux density. Using (5.1) -(5.3) we obtain 



(5.7) 



r=r+ 



fJZ n-i /n~l n-l /m-2„ 



for the temperature, entropy density, and horizon fiux density. 
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5.2 Integrated First law of thermodynamics 

Before embarking on our numerical calculations, in this section we obtain relationships between the free 
energy density and the energy density derived in section 4 and the above thermodynamic quantities. We 
then use this to obtain analytic predictions when the marginally relevant mode vanishes, i.e. A — t- 0. 

First, by using the asymptotic solutions (3.1) - (3.3) we construct an r-independent RG-invariant 
quantity 



K -- 

2 



-\/ f{r)p{r) 2 ^ — g(r) + m(r) + A;(r)x(r)^ , 



n-l 

(r) 2 



4x + 4(n — l)zk{r) — n{n — l)q{r) — x{r) + 2{n — l)q{r)k{r)x{r) I , (5.9) 



4(n — l)q{r 

where plugging (3.6)-(3.8) into the above gives 

V^^(l-2z)V (l-2z)(z-l)2'^y ^ ^ 

Near the horizon, this RG-invariant quantity can be calculated by using (5.1) -(5.3) and expressed in 
terms of T and s by applying (5.8). We find 

K=^,p^ =Ts^. (5.11) 
Next, from the free energy density (4.8) and the energy density (4.19) we obtain the following relation 

which is more simply expressed as 

^n— 1 

£ = J^+ ^i^. (5.13) 

This relation is easily checked using (4.24), (4.25) and (5.11). 
Finally, combining (5.11) with (5.13) gives 

F = £-Ts. (5.14) 

which is the integrated form of the first law for these black holes. We will use this to check the accuracy 
of our numerical results in section 5.4. 

Considering the limit A = 0, since anisotropic scale invariance still holds, from the Ward identity we 
expect that the pressure is equal to the energy [12]. From (4.20), we have 

Jo = -^o (A = 0), (5.15) 

and in conjunction with (5.14), we obtain an analytic prediction for when the marginally relevant modes 
are not excited 

Fo = -£q = -\ts^ (A = 0), (5.16) 

which will be used for a consistency check on our numerical results in section 5.4. Note that the relation 
T = —8 also holds for A ~ case when /3 = 0, as is easily seen from equations (4.24) and (4.25). 
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5.3 Integrating towards the Lifshitz Boundary 

To investigate Lifshitz spacetime in the UV-region, T ^ A^, we re-expand the asymptotic solutions into 
the high energy regime by using (3.5), and set the arbitrary scale ^ ~ r^^, thereby making the marginally 
relevant modes near-Lifshitz. Upon carrying this out, the deformed Lifshitz spacetime is described in the 
high temperature regime with fixed A ~ 0. This spacetime approaches a pure-Lifshitz by supplying heat at 
much greater temperatures or higher energy scales. This process of expansion yields a constant log(Ar+) 
which must be fixed in the asymptotic regions. Near the horizon the expanded black hole solution has 
three parameters /o, po, and ho where /o and po will be determined, depending on which asymptotic 
spacetime is imposed. 

In this section, we numerically connect the expanded near-horizon black hole solution with the near- 
Lifshitz asymptotic solution. By matching two solutions, we determine the constant log(Ar_|_) which 
corresponds to ho in section 5.3.1, and fo and po in section 5.3.2. Then with these values we compute 
the thermodynamic quantities, varying the value of ho- We compute the entropy density as a function 
of log(A^/T) in section 5.3.2, and the free energy density and energy density versus log (r/r+) in section 
5.3.3. In section 5.4, we compute the free energy density and energy density as functions of log(A^/T) 
and find a suitable fitting curve, obtaining a prediction on the sub-leading order of the free energy and 
energy density as a function of log(A^/T). We also discuss how renormalization group flow is described 
in our context in section 5.3.1. In our numerical work, we use r/r+ as our radial variable, and unitless 
quantities such as T/Ts, £/Ts, and A^/T are considered. 

5.3.1 Matching A 

For extracting A, we start with the asymptotic solutions k{r), q{r), and x{r) in (3.1) - (3.3). As the 
UV-deformation is applicable to Lifshitz spacetime in the high energy limit A^/T — )• 0, we expand the 
solutions as in (3.6) - (3.8), bringing in the arbitrary scale /x >> A eventually setting fj, — )• r^^ . Then, 
near the asymptotic boundary, we have only to fix the variable log(Ar+). Near the horizon, the expanded 
black hole solutions (5.1)-(5.3) transfer to functions of k{r), q{r), and x{r) by (2.15) and (2.16), and in 
this process fo and po drop out and ho only remains. Now, having a proper value of ho, we numerically 
integrate from the near horizon to the boundary satisfying the equations of motion (2.17), matching the 
numerical solutions to the asymptotic expectation by adjusting the value of log(Ar4-) in the middle region. 
In our numerical works, the integration starts at log(r/r_|_) ~ —0.015 and ends at log(r/r+) ~ —10^. Our 
numerical results for the k function are shown in figure 1, where the red dashed line is a fit for the 
asymptotic expectation and the blue solid line corresponds to the numerical results. We note that the 
agreement between our numerical results and the asymptotic expectations is very strong. 

The important outcome of this procedure is that a maximum value of ho is obtained. In other words, 
if we keep increasing the value of ho then beyond a certain point we are not able to find a matching 
condition connecting our numerical result to the asymptotic expectation. Physically this means that at 
a large value of fiux, the metric functions grow exponentially as the boundary is approached, so they do 
not ever reach the boundary. This means that under the condition of large fiux and high temperature 
A^/T — )• 0, the spacetime having the black hole (5.1)-(5.3) is no longer deformed, but rather is asymptotic 
to pure Lifshitz spacetime. We also find that the maximum values of ho linearly increased according to 
the critical exponent, z (or the spatial dimension) of the spacetime, n. We present this behaviour in 
Figure 2, and explicitly denote the maximum values of ho in Table 1. 

On the other hand, a minimum value of the fiux also exists and its value becomes zero at the horizon. 
In this case, the massive vector field disappears and the spacetime is described by an asymptotically AdS 
black hole. Examining the situation for ho between h^ax and hmin is non-trivial. However for a small 
amount of fiux we can verify that in the zero temperature limit A^/T — )• oo the renomalization group fiow 
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log 



(a) z = 2 and ho = 0.97128, log(Ar+) = -11141.7 



(b) z = Z and ho = 1.63428, log(Ar+) = -4188.2 



log 



(c) z = A and ho = 2.28218, log(Ar+) = -2709.4 



0.894403 
0.894403 
0.894403 
0.894403 
0.894403 
0.894402 
0.894402 
0.894402 



(d) 2 = 5 and K = 2.92548, log(Ar+) = -2270.6 




(e) z = 6 and ho = 3.56678, log(Ar+) = -2726 (f) 2 = 7 and ho = 4.20680, log(Ar+) = -908.7 

Figure 1: Extracting log(Ar+) 
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Figure 2: /imax versus z 



is recovered [4] . This obviously indicates the existence of an RG flow with the marginally relevant mode, 
and also implies that tuning the horizon flux via ho interpolates between the zero temperature RG flow 
[4] and asymptotic Lifshitz black holes [7, 12, 13, 14, 15, 16]. 

5.3.2 Matching /o, and pQ 

Since, given /iq, log(Ar+) is fixed, /o and po arising in (5.1) - (5.3) can be determined. Repeating the 
previous process with the functions of /, p, and h, the asymptotic solutions (3.9) - (3.10) are in the high 
energy regime obtained via the expansion of (3.5) and the rescaling of (3.13), becoming the expressions 
(3.11)-(3.12). We can neglect the a and /? involved in the exponent terms because these contribute much 





hmax 


z = 2 


0.9713 


z = 3 


1.6343 


z = A 


2.2822 


z = 5 


2.9255 


z = Q 


3.5668 


z = l 


4.2070 


z = 8 


4.8465 


z = 9 


5.4856 


z = 10 


6.1244 


z = ll 


6.7629 


z = 12 


7.4012 


z = 13 


8.0394 







Table 1: maximum value of /iq 
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(a) For z = 2 and ho — 0.962, which corresponds to log(Ar-+) = —106.8, the red dashed hne is matched to the blue soHd line 
at forX = 29.04 (left) and at por^ = 1.8428 (right). 




(b) For 2 = 3 and ho = 1.628, which corresponds to log(Ar+) = —72.6, the red dashed line is matched to the blue solid line 
at fori = 44.266 (left) and at por+ = 1.7579 (right). 




(c) For z = 4 and ho — 2.2798, which corresponds to log(Ar+) — —133.4, the red dashed line is matched to the blue solid 
hne at fori = 71.40 (left) and at porl = 1.5781 (right). 



Figure 3: Extracting /o and po for z = 2,3, and 4 
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(a) For z = 5 and ho = 2.9238, which corresponds to log(Ar+) = —155.6, the red dashed hne is matched to the blue solid 
line at for^ = 109.953 (left) and at por+ = 1.45079 (right). 




(b) For z = 6 and ho — 3.5655, which corresponds to log(Ar+) — —183, the red dashed line is matched to the blue solid line 
at /or+ = 157.777 (left) and at por+ = 1.36751 (right). 




(c) For z = 7 and ho = 4.206, which corresponds to log(Ar+) — —220.2, the red dashed line is matched to the blue solid line 
at for+ = 214.442 (left) and at por+ = 1.3098 (right). 



Figure 4: Extracting /o and po for z = 5,6, and 7 
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less than the logarithmic terms as the bomidary is approached (technically these unknowns are hard 
to fix). Applying the same numerical technique, we start to integrate from the near horizon towards 
the boundary governed by the equation of motion (2.14). Over a finite range, provided values of ho 
and log(Ar+), the numerical solutions and asymptotic expectations (3.11)-(3.12) are matched by finding 
values of /or_|_^^ and po'^+^- Our numerical calculations are carried out over a range of 4 to 9 dimensional 
spacetime. The results are shown in figures 3 and 4 where the red dashed line is a plot of the asymptotic 
expectation and the blue solid line is the numerical result. 

Using the fixed constants A, /o, and po we next explore physical quantities of interest. Computing 
the entropy density, derived in (5.8) and determined by /or+^^ and por+'^, we plot s/T as a function of 
log(A^/r) in figure 5. Here, the value of log(A'^/r), equivalent to log(Ar+), is obtained from the relation 
between the thermodynamic variables in (5.8), which is 



We explicitly consider dimensionalities ranging from 4 to 9. 
5.3.3 Energy Density, and Free Energy Density 

Now we move to the free energy density (4.8) and energy density (4.19). Putting the numerical results 
of k, q, and x functions into (4.8) by using (2.15) and (2.16), and applying the counterterm (4.21) - 
(4.23), the free energy density and the energy density are respectively depicted in figures 6 and 7. We 
find that the free energy density and the energy density have a flat region over some finite range of 
log(r/r_|_), yield a stable constant value for these quantities. In these figures we also illustrate oscillating 
and divergent behaviours as the boundary is approached. Recall that initially the physical quantities 
having marginally relevant modes diverge at the boundary, necessitating the addition of counterterms 
to render them finite. The counterterms should be expanded in an infinite series in log(Ar). However 
in practice the counterterms (4.21) - (4.23) are truncated at a finite order, while our numerical results 
include higher orders than ones we considered in the analytic calculation. This limitation is responsible 
for the unstable behavior near the boundary that appears in figure 6 and 7. 

As the spacetime dimension increases we find that the flat region gets narrower and the divergence 
behaviour starts more quickly from the horizon. The main reason for this is due to the \/7p^^^ factor 
commonly appearing in the free energy density (4.8) and the energy density (4.19); both ^/f and p^^'^ 
vary as 1/r^ (as shown in (3.12)) as the boundary r — )• is approached, whereas the quantities in the 
parentheses of (4.8) and (4.19) contain no divergent terms. 

5.4 Exploring the dependence of £ and T on logA^/T 

In this section, we numerically compute the free energy density per Ts (J-'/Ts), and the energy density 
per Ts {£/Ts), as functions of log(A^/T), and investigate their behaviour in terms of T, assuming that 
A is very small and flxed. Figure 8 and 9 present our results for these quantities and for J^/£. We also 
include the fltting function depicted as a solid line. As expected from analytic considerations for the 
leading order terms of J-'/Ts and £/Ts from (5.14) for A = 0, we recover from our numerical results 
the same value for the leading order terms. In addition, as the marginally relevant modes generated by 
A ~ are numerically evaluated, their effect is also shown in Figure 8 and 9. Their behaviour predicts the 
sub-leading terms expressed as functions of log(A^/r). For each z, the fltting functions of J-/Ts, £/Ts, 
and J^/£ are presented in Table 2. 




(5.17) 
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(a) Dots are calculated under that z = 2 and ho runs from 
0.9713 to 0.9705 in increments of 0.00004 and are joined 
by straight lines. 



11.90 - 

T I* 
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(c) Dots are calculated under that z = A and ha runs from 
2.2822 to 2.2814 in increments of 0.00004 and are joined 
by straight lines. 




-20000 -15000 -10000 -5000 




(e) Dots are calculated under that 2 = 6 and ho runs from 
3.5668 to 3.566 in increments of 0.00004 and are joined by 
straight lines. 




-14000 -12000 -10000 -8000 -6000 -4000 -2000 




(b) Dots are calculated under that z — 3 and ho runs from 
1.6343 to 1.6335 in increments of 0.00004 and are joined 
by straight lines. 
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(d) Dots are calculated under that z — 5 and ho runs from 
2.9255 to 2.9247 in increments of 0.00004 and are joined 
by straight lines. 



6.955 



6.940 




-25000 -20000 -15000 -10000 -5000 




(f ) Dots are calculated under that z — 7 and ho runs from 
4.2070 to 4.2062 in increments of 0.00004 and are joined 
by straight lines. 



Figure 5: Entropy density per unit temperature versus log(A^/T) 

19 



(a) When z = 2 and ho = 0.97102, corresponding to log(A^/r) ~ —5270, the left is the numerical result of free energy 
density F over Ts and the right is the numerical result of energy density £ over Ts, as a function of log(r/r+). The quantity 
is well defined in the intermediate region. Red dashed line is reading-off the constant value of T jTs in the intermediate 
regime. 




(b) When 2 = 3 and fto = 1.63422, corresponding to log(A^/r) ~ —8000, the left is the numerical result of free energy 
density T over Ts and the right is the numerical result of energy density E over Ts, as a function of log(r/r+). The quantity 
is well defined in the intermediate region. Red dashed line is reading-off the constant value of T (Ts in the intermediate 
regime. 




(c) When z — 1 and /lo = 2.28212, corresponding to log(A'^/T) ~ —7000, the left is the numerical result of free energy 
density T over Ts and the right is the numerical result of energy density E over Ts, as a function of log(r/r+). The quantity 
is well defined in the intermediate region. Red dashed line is reading-off the constant value of T jTs in the intermediate 
regime. 



Figure 6: F /Ts and £ /Ts versus log(^) for z = 2, 3, and 4 



(a) When 2 = 5 and ho — 2.92526, corresponding to log(A^/T) ~ —4100, the left is the numerical result of free energy 
density T over Ts and the right is the numerical result of energy density £ over Ts, as a function of log(r/r+). The quantity 
is well defined in the intermediate region. Red dashed line is reading-off the constant value of T (Ts in the intermediate 
regime. 




(b) When 2 = 6 and — 3.56670, corresponding to log(A^/r) ~ —8700, the left is the numerical result of free energy 
density T over Ts and the right is the numerical result of energy density E over Ts, as a function of log(r/r+). The quantity 
is well defined in the intermediate region. Red dashed line is reading-off the constant value of J- /Ts in the intermediate 
regime. 




(c) When z — 7 and ho — 4.20694, corresponding to log(A^/T) ~ —14000, the left is the numerical result of free energy 
density F over Ts and the right is the numerical result of energy density £ over Ts, as a function of log(r/r+). The quantity 
is well defined in the intermediate region. Red dashed line is reading-off the constant value of T (Ts in the intermediate 
regime. 



Figure 7: F /Ts and £ /Ts versus log(^) for z = 5, 6, and 7 



In addition, the marginally relevant mode should be consistent with the first law of black hole ther- 
modynamics, which is —J-/Ts + £ /Ts — 1 = 0. We plot this in figure 10 as a check on the accuracy on 
our numerical results. We find that our numerical errors are found between the order of 10~^ and 10~^. 





free energy density over Ts 
(ts =) 


energy density over Ts 

(ts =) 


free energy density over 
energy density =) 


z = 2 


1 1 1 

9 Irio- /T ~^ 
^ il /J. 


1 1 1 

9 Ino- /T 
^ l'-J& *- / -'- 


1 4 1 

log ii / J 


z = 3 


1 0.76 1 
i logA^/T 


1 0.76 1 

2 logA'^/T 


1 3 1 
logA'^/T 


z = 4 


1 0.67 1 

2 logA4/T 


1 0.67 1 

2 logA-i/T 


1 2.7 1 

^ logA4/T 


z = 5 


1 0.63 1 

2 logA5/T 


1 0.63 1 

2 logA5/T 


1 2.50 1 
^ logA5/T 


z = 6 


1 0.60 1 

2 logAB/T 


1 0.60 1 

2 logA6/T 


1 2.40 1 
^ logA6/T ^ 


z = 7 


1 0.58 1 

2 logA'^/T 


1 0.58 1 

2 logA7/T 


1 2.33 1 
^ logA^/T 











Table 2: fitting functions for Ts' i 



6. Summary and Discussion 

We have carried out a study of the deformations of Lifshitz spacetime in higher-dimensions, where the 
deformation is produced by the marginally relevant modes due to A ~ 0, which is the dynamically 
generated scale. Our main objective was to describe the renormalization group flow of the marginally 
relevant operators and to investigate the thermodynamic behaviour at flnite temperatures in terms of 
the dimensionality of the spacetime. We began with the assumption [4] that the deformation is small 
in the UV-region and leads to renormalization group flow into IR-region, where the Lifshitz spacetime 
is described in the high temperature limit A^/T — t- and the AdS spacetime is formed in the zero 
temperature limit A^/T — )• oo. 

We derived the equation of motion for Einstein gravity with the massive vector fleld, and set up the 
ansatz for the metric and the vector potential for both Lifshitz spacetime and AdS spacetime in section 2. 
In section 3, we obtained asymptotic solutions in terms of functions k, q, and x and also in functions of /, 
p, and h (used separately for later numerical work), and derived the free energy density and the energy 
density £ at the spacetime boundary in section 4. Since the marginally relevant modes make the action 
and the physical quantities diverge, we performed holographic renormalization by adding counterterms, 
constructed from the covariant quantity (fc^ — (^z — l)/z). The coefficients Cj in the counterterm expansion 
were obtained as a series in log(Ar+). We also obtained near-horizon expansions of (planar) black hole 
solutions in section 5.1, and proved that our analytical results are consistent with the first law of the black 
hole thermodynamics as A = 0, which means the marginally relevant modes are turned off. Numerically 
we computed in space-time dimensions 4 to 9 the entropy density s/T, the free energy density T /Ts, and 
the energy density E /Ts as functions of log(A^/T) in near-Lifshitz spacetime. 
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Figure 8: Plots of T jTs , E jTs and T jE versus log(A^/T) for 2 = 2, 3, and 4. 
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(a) Plots of T/Ts, £/Ts, and J-/£, as a function of log(A^/r) for z — 5. Dots are numerical results running ho from 2.9255 
to 2.9249, which corresponds to log(A^/T) from about —14000 to —2000. Solid line is fitting equations in table 1. 
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(c) Plots of T/Ts, £/Ts, and T/£, as a function of log(A^/r) for z = 7. Dots are numerical results running ho from 4.2070 
to 4.2064, which corresponds to log(A^/T) from about —30000 to —2200. Solid line is fitting equations in table 1. 



Figure 9: Plots of T jTs , E jTs and T jE versus log(A^/T) for z = 5, 6, and 7. 
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(a) Plots of -F/Ts + £/Ts - 1 for 2 = 2 and feo from 0.9713 
to 0.9707. 



-0.0000361 
-0.00003611 
-0.00003612 
- - 1 -0.00003613 
-0.00003614 
-0.00003615 
-0.00003616 



-14000 -12000 -10000 -8000 
I A? 



-6000 -4000 -2000 



(b) Plots of -F/Ts + £/Ts - 1 for z = 3 and ho from 1.6343 
to 1.6337. 
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(c) Plots of -F/Ts + e/Ts - 1 for 2 = 4 and ho from 2.2822 
to 2.2816. 



F E 

Ts Ts 



-0.00010860 
-0.00010865 
-0.00010870 
-0.00010875 
-0.00010880 
-0.00010885 
-0.00010890 



-12000 -10000 -8000 -6000 -4000 -2000 



(d) Plots of -F/Ts + - 1 for 2 = 5 and ho from 2.9255 
to 2.9249. 
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(e) Plots of -F/Ts + f/Ts - 1 for 2 = 6 and ho from 3.5668 
to 3.5662. 
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(f) Plots of -F/Ts + f/Ts - 1 for 2 = 7 and ho from 4.2070 
to 4.2064. 



Figure 10: Entropy density over T versus log(A^/T) 
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Our results indicate that the basic physics of Lifshitz/QCT duahty [6] is vahd in higher dimensions. 
Regardless of dimensionality, with a small flux ho, renormalization group flow under the marginally 
relevant operators for UV-Lifshitz to IR-AdS was obtained in the zero temperature limit A^/T — )• oo, 
commensurate with the 4-dimensional case [4]. This implies that we can expect RG flow to exist for each 
horizon flux, ranging from the zero temperature limit to asymptotically Lifshitz black hole spacetimes. The 
thermodynamic quantities, s/T, T/Ts and £/Ts were calculated as functions of log(A^/T) just below 
the maximum value of Hq (which describes a very slightly deformed Lifshitz space-time) and showed 
graphically similar behaviour. Also, as analytically expected, the leading-order terms from T jTs and 
£/Ts from the first law of black hole thermodynamics for A = can be numerically obtained. We 
also found the sub-leading terms in terms of log(A^/r) when A ~ 0. This illustrates how physical 
quantities change as functions of log{A^)/T upon approaching the critical point from a given phase. From 
a holographic duality perspective we can characterize the strongly coupled physics having the marginally 
relevant mode. 

Since we considered the case for which the energy- momentum tensor of the gravitational field and the 
operator of the vector field become marginally relevant, we set z = n — 1. We consequently found that 
higher dimensional Lifshitz black holes can hold more flux, since the maximum value of ho increased with 
increasing n. The sub-leading term caused by the marginally relevant modes on the J-'/Ts and £/Ts has 
a weaker dependence on temperature T as dimensionality increases, as shown in table 2. 

There are several possible directions for further research. One involves studying the relationship 
between the RG flows described here and previously noted stability issues in asymptotically Lifshitz 
space times [17]. Another involves consideration of marginally relevant modes in Lifshitz spacetimes with 
higher curvature corrections. The action with a Gauss-Bonnet term is now in progress [18]. 
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